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Detailed BT Data for ICAPS 
(plots reproduced from reference 2) 
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SIMAS Historical Profiles (reference 2) 
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BT Data for SIMAS 
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APPENDIX A 


Section A6 
Sound Speed Profiles (reference 2) 
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Historical SSP Plots (reference 2) 


Solid lines identify ICAPS profiles 
Dotted lines identify SIMAS profiles 
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APPENDIX D 


ICAPS and SIMAS Sound Speed 


Profile Generation Methodologies 


MEMORANDUM 


TO: LCDR AI Galus 

FROM: J. Locklin and B. W. Scaife 

DATE: 21 April 80 | 

SUBJECT: ICAPS and SIMAS Sound Srecd Profile Gererz tion Methodolo- 
gies 


ENCLOSURES: (1) Features of ICAPS and SIMAS SSP GENERATORS 
(2) ICAPS Logic Flow 
(3) SIMAS Logic Flow 
(4) Problem Areas in ICAPS and SIMAS Codes 
(5) ICAPS Temperature Syithesizing Algorithra 
‘6) Test Resu'ts 


REFERENCE: (1) "The ICAPS Water Mass History File," by Alvan Fisher, 
Jr., NOO RP-~19, May 1978. 


(2) "Description of ICAPS Environmental Data Structure" by 
John Lever, NAVOCEANO TN 3700-82-79, March 1979 


I. The purpose of this task was to provide tne capability to examine the 
parameter (temperature, salinity, depth, sound speed) and coded algorithm 
depenuencies of the SIMAS and ICAPS proiile generation processes. This 
effort is directed toward determining an optimal choice of parameters to 
be included in the APF Data Base which is to support both ICAPS and 
SIMAS as they are to be configured in the Comprehensive APP systems. 


2. Enclosures (1), (2), and (3) document the respective algorithms and were 
prepared from FORTRAN source code listings. ICAPS listings were 
provided by Mr. Doug Gordon, NAVOCEANO, for both the Univac 1108 and 
NOVA 800 software versions. SIMAS listings were provided by Mr. Burley 
Brunson, NORDA CODE 320, for tue PDP 11/4u version instalied at 
NORDA. In the process of examining and documenting the source code, 
several problems were discovered in each of the systems. These problems 
are discussed in enclosures (4) and (5). 


3. In order to study the performance of each algorithm, a version of each was 
coded from the source listings and installed on the CDC 6€00 for test 
purposes. The basic algorithms to treat historical and 37 projsiie data, and 
to prepare resultant SSP's remain intact and reflect the same methods used 

. in ICAPS and SIMAS codes currently in operation on their respective 
computers. 


| 
} 


4, A preliminary set of tests were conducted to verify thc CDC versions, look 
at the criticality of near-surface vertical salinity variaiion, and compare 
the merge processes. E. hiashimoto, NORDA code 321, provided hydro- 
graphic data at two stations. The results of these tests are contained in 
enclosure (6). The test cases illustrate the applicability of the ICAPS and 
SIMAS algorithms. 
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i ic 


Se 


ENCLOSURE | 


FEATURES OF iCAPS AND SIMAS SSP GENERATORS 


Enclosure (1) to 
(OSI memo dtd a/21/8y 


Jocklin to Galus 


METHODOLOGICAI. FEATURES OF SSP GENERATOR 


Multip!e water mass (up to 5) capability 
for each geographic area. 


Histurical data base contains for each 
water mass, temperature and salinity 
profiles (2H, THis SH;), izl, NH «45. 


Extends the Bathy thermography (BY) 
profile, (ZB., TB.), j-J, NBT; vith 
modified higtorickl temperatures jor 
corresponding water mass at depths be- 
low the last BT depth, and intetpolate 
salinity values from historical profile. 


Insert bathymetry and truncate SSP at 
bottom depth. 


Correct salinity values for a stable en- 
vironment. 


Utilize Wilson's equation to calculate 
SSP. 


Calculate the layer depth. 


| Calculate the layer depth. 


fe lo Univac 
(1) NBi fi Raga 


A geographic area is represented by 
one water mass. 


Historical data base contains for 
each water mass a sound speed pro- 
file (Zi, VH;), NH €27, and a rep- 
resentative near-suriace _ salinity 
value, S. 


Utilize Leroy's equation to calculate 
the SSP for BT profile with histori- 
cal constant salinity. The BT pro- 
file may consist of up to 24 points. 


Shift the BT SSP by the lin2arly 
interpolated velocity difference be- 
tween historical SSP and BT SSP at 
1000 ft. Next, extend below the BT 
SSP with unmodified historical SSP 
points. 


Insert bathymetry and truncate SSP 
at bottom depth. 


Adjust fathometer bathymetry and 
interpolate sound speed at adjusted 
bottom depth. 


ENCLOSURE 2 
ICAPS LOGIC FLOW 


Enclosure (2) to 
ODSI memo dtd 4/21/80 
Locklin to Galus 


ICAPS Sound Speed Profile (SSP) Generation (UNIVAC Version) 


Input identifying data 
(date, bathymetry, 
geographic position 
and a water mass 
selection number, etc.) 


Ts bathytnarmegraph 
(BT) data to be input? 


No Yes 


Input BT data (denth, 
temperature) NBT pairs 
of profile points (4), 
Corvert to metric if 
indicated, 


Retrieve specified 
historical orofiie 
(depth, temperature, 
salinity) for this 
water mass. If no 
bathymetry was input, 
also retrieve histori- 
cal bathymetry. 


Is this a shaliow 
water case? (Maximum 
BT depth less than 
200 meters). 


Define the estimated 
profile (depth, 
temperature, salinity) 
to be the historical 
profile, and terminate 
the prorile at the 
linearly interpolated 
bathymetry. 


Define the estimated 
profile (depth, tem- 
perature, salinity) 
so that the (depth, 
temperature) profile 
is the BT profile and 
the sal‘nity values 
are corstant (35PPT). 
Set the bathymetry to 
be the maximum BT 


ae | | depth. 


16 Ne 


(1) NBT € {35 Nova D-6 


eee eS ER = 


._—- 


Select historical profile (depth, temperature, salinity) for this water mass 
~~ NH number of triple points and retrieve the merge weighting factor (MFACT). 


Retrieve up to five 
historical prufiles 
with 200m temperature 
tolerances and 

200-306m gradient 
tolerances. 


Linearly interpolate 
the 200m and 300m BT 


temperature yalues and 
detine the 200-309%m 
BT cradient. 


Is there more than one 
historical profile? 


Select first histori- 
cal profile. 


Has an historical pro- 
file been specified? 


Setect specified 
historical profile. 


Select an historical 
profile using 


jbo Jen ADEES: 


[ts bathymetry needed? 


Yes 


| Retrieve historical 
| bathymetry. 


me 


Is any BT profile 


depth 2 the 
bathymetry? 


Is the maximum BT 
depth greater than 50 
meters above the 
maximum historical 
profile depth? 


No | Yes 


(Find index I to ter- 


minate BT profile at 
the shallowest BT 
depth = bathymetry. 
Set NBT=I. 


Define estimated pro- 
file (depth, tempera- 
ture) to be truncated 
BT profile. 


Define the estimated 
profile (depth, tem- 
perature) to be the 

BT profile. 


—_—_— 


The estimated profile 
J has N=NBT pairs of 
points (2M;, TM; ) 


Select the index, k, 
for the shal lowest 
depth, Z,, in the 
historical profile 

where the maximum 

BT depth is «€ 

Z,-50 meters. 


Initialize the first 
NBT number of genera- 
ted (depth, tempera- 
ture) profile points 
to be the BT profile. 


Complete the (depth, temperature) profile being generated by augmenting the 
initialized NBT number of points with the historical depths Zj for i=k,...,NH and 
modify the associated historical temperatures, T 


eee mmm 


(TT, 4) (ZMy-Z, 
Calculate aT=TM- 4+ — spo k-1 
i 


X=MFACT/1000. 


Loop over historical profile for i=k,...NH 


oe 


Calculate N =N+1 


ZMy=2 5 


at =( ar)x (24725 ~1)/100. 
TMy=T 5 + OT 


Is the depth, iM), 
less than the 
bathymetry? 


Is Weopee comp leted? 


The estimated profile 
now has N pairs of 
points 


Calculate the estimated salinity values, SM; at the depths 
ZM., i21,...,N. 
see s . 


——  ——_—_ — — — mea esa eee ele 


me ae eee ie 


Is the estimated depth 
ZM., less than the 
deepest historical 
depth? 


Linear ly interpolate 
salinity from the 
historical profile. 


Set salinity constant 
as the last interpo- 
lated value. 


Is loop compieted? 
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Update the estimated profile so that the deepest depth is 
the bathymetry. 


Is batnymetry greater 
than ZMy? 


Linearly interpolate 
temperature and 
Salinity, id set 
ZMy=bathymeury. 


Linearly extrapolate 
temperature and 

salinity, and add a 
point to the profile 
at ZM.4,=sathymetry. 


Set N=N+1, 


Correct the estimated salinity values for .tability. 


Calculate the densities by the density 
funccion OF (TM; ,SM;), for i=1,..., NBT. 
j 


Loop ove over jel =), NET-4, | 


_ 407-0.01 
.77184(1.69)1 69)19°/ £¥ 


_ over i= Phaexe aN. 


= 


Calculate sound speed 
at depth ZM, via 
Wilson's equation. 


Is loop a 


~ 


Calculate the Layer 
depth. 


EXIT 


ICAPS sound speed 
profile generation 
comp leted! 


ENCLOSURE 3 
SIMAS LOGIC FLOW 


Enclosure (3) to 
ODSI memo dtd 4/21/80 
Locklin to Galus 


SIMAS Sound Speed Profile (SSP) Generation 


Input identifying data. 


Is SSP generation 
method to consider 
bathythermocraph (BT) 
data or previously 


ated nrofile? 
Yes No 
ENVIRN FORCST 


eens 3 


Inoue latest SSP Input MGS providence 
generated. chart depth or 


fathometer reading, 
and wind speed. 


Input MGS providence, 
chart depth or 
fathometer iveading, 
and wind spead. 


. STcerses Is there manual SSP 
used? 
Retrieve historical - Input SSP data points - 
SSP and representative (depth, sound speed) - 
sclinity value: (Z;,C;), i=1,N « 50, 


((ZH,,VH,), =1,NH),S. 


Retrieve historical 
SSP and representative 


salinity value: Set bath 
rt ymetry to Z(N). 
((2H;,VH;), j=1,NH),S. 


Set SSP being genera- 
ted to be the histori- 


oS — cal SSP. 
Insert the historical, Let N =NH, 
linearly interpola‘ed, Ss ee ce 
sound speed point at | zs ZH; | J=1.N 
1000 ft. into the | C.2VH 
historical protile. es aa | 


7 
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Is BT data to be used? 


[ input BT data (depth, Set SSP being genera- 
tenperatur2) pairs: ted to be the histori- 
(ZB., TB;), i=l,NBSX4. cal profile. 


Are any points to be 
edited? 


Input d-ta to add, de- 
lete or change any 
point. 


METRIC input units. 


ae Gre 


Convert BT points to 
English units. 


Via LEROY's equation, 
calculate sound speeds 
(VH,) usirg BT data 
and historical salini- 


ty. VH = p (2B; ,TB;,S) 


Calculate the layer 
depth. 


Insert the souna speed 
point at i000 ft. ino 
the BT SSP. 


feet? 


Soeeeee | 


| insert the lingarly 


interpolated point at 
1500 ft. into the 
nistorical SSP. 


Does user want to in- 
put another BT profile? 


Are the sound speed BT 
and historical alues 
within 9 tc./sec. of 
each other at 1000 


Revise the BT SSP to 
be a modified histori- 
cal SSP down to 1500 
ft., so that the layer 
deptn is the same as 


the initial BT Sr 
| layer depth. 


Revise the BT SSP by 
shifting the initial 
BT SSP so that at 1000 
ft. the 37 sound 

speed is the same as 
the historical value. 


Initialize the SSP 
being generated to 
consist of ihe revised 
BT SSP points. 


' Extend the initialized 
| new SSP with those 

historical SSP points 
that are below the 
maximum BT depth. 


D-17 


| oO ss 
Insert the linearly 
interpolated bathym- 


etry point into the 
generated SSP, 


eee aes 


Truncate the genera- 
ted SSP at bathymetry 
[i.e. (Z., C;) i=1,N 


& Z,=bathymetry] . 


Lo:p for i=2,N. 


cates —_— — — — — —, ~— 


> 


Loop done? 


No Yes 
| Set C.=C-+0.01. | 
| 
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Perform bathymetry correction: 


4 isn 
Set SUM = Tl, PS (2,-2.1) (C;+C,_3) 
Zn 
ZB = agg SUM 


at corrected depth cB, linearly 
interpolate sound speed and replace 
(Zy, Cy) by this new corrected point. 


Calculate the layer 
depth. 
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ENCLOSURE 4 


PROBLEM AREAS IN ICAPS AiND SIMAS CODES 


During the implementation process to establish ICAPS and SIMAS on the CDC 
6600 computer for comparative analysis purposes, problems were revealea. 
These problems are assembled for each program under three categories: 1. 
program errors, 2. program cautions, and 3. cominents. 


Enclosure (4) to 
ODSI memo dtd 4/21/80 
Locklin to Galus 


1.1.1.2 


1.1.1.3 


1.1.1.4 


L.L1.5 


ICAPS 
PROGRAM ERRORS 
Univac Version 


The synthesized profile contains salinities that are initialized by 
interpolating values from the historical profile, and wh2n a depth of 
the synthesized profile is equal to or exceeds the maximum depth of 
the historical profile (a most frequent occurrence) the salinity is sci to 
be the historical salinity for the previous depth. This can be corrected 
by changing the inequality ".LT." to ".LE." in the IF statement 
following FORTRAN statement number 50 in subroutine MERGE. This 
code error will have negligable impact for the deeper profiles. For 
shallower profiles the irnpact will depend upon the effect on velocity 
by the error in salinity. 


When the user wishes to use only the historical data (no Bi daia: 
NDP=0 and BOTZ 1.0) and metric is not specified (MOE £ "M"), the 
program incorrectly multiplies the historical bathymetry in meters by 
the feet-to-meters conversion factor, FTMT. This can be corrected by 
inserting after format statement number 9001 in the main program 

IF (MOE .NE. 1 HM) BOTZ = FTMT*BOTZ 


and deleting the statement immediately after FORTRAN statement 
number 10, and replacing FORTRAN statement number !60 with 
100 CONTINUE. 


An error message is not utilized when appropriate and an incorrect 
message is printed. Change the line just preceeding format statement 
number 1045 to read 

WRITE (6, 1045) IBTYP2. 


In the main program, there is an "MP" variable occurring just after 
format statement number 1013 that is not initialized. Apparently, this 
was to have been defined by "NP" in the preceeding call to subroutine 
COMPUX. This is cerrected by changing "NP" to "MP," the Jast 
argument in the argument list for the "CALL COMPUX ..." statement, 
two lines after FORTRAN statement number 360. 


Subroutine MERGE calls subroutine SGMTST to perform the water 
column stability algorithm trom the bottom of the merged profile. 
Though this is in contrast with the NOVA version that starts from the 
bottom of the BT, it is more reasonable since the algorithm is applied 
just below the BT where merging occurs (refer to paragraph 1.1.2.2). 
Subroutine SGMTST has an incorrect DIMENSION statement allocating 
insufficient memory, particularly for the density array. The current 
dimension for 30 elements is not consistent with the calling routine. 
The arrays should be dimensioned for 50 elements each, since the 
calling routine, MERGE, identifies the arrays as for the merged 
profile. This is corrected by changing the dimension statement in 
subroutine SGMTST to read 
DIMENSION T(50), S(50), SIGMA(50) 
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1.1.1.6 


L.1.1.6.1 


1.1.1.6,2 


1.2 
1.2.1 
1.2.1.1 


1.2.1.2 


1.3 


Under normal operating system usage, the density and stability calcu- 
lations will be incorrect. Mixed expressions contain integers raised to 
negative integer powers which result in zero values. 


Subroutine SIGMAT has an equation "B =..." that contains the integer 
10 raised io a negative integer power. 


Subroutine SGMTST has the equation "DS =..." that contains two 
integer 10's raised to negative powers. 


Nova 800 version 
Same error as described in paragraph 1.1.1.1 above. 


The progratn should perfoiin the wate: column stability algorithm 
from the bottom of the merged profile as the Univac version does 
(refer to paragraph 1.1.1.5) -- rather than from the bottom of the KT. 
Thus, the ejfects on density (refer to paragraph 1.3.4.2) due to 
temperature revisions by merging just below the deepest BT depth 
would be accounted for. his can be corrected by changing in 
subroutine MERGE, the FORTRAN statement just after statement 
number 100 to: 
N=NOPTM- 1 


and two lines later, change tu read 
K=NOPTM-J. 


PROGRAM CAUTIONS 
ICAPS, Univac and Nova versions 


When a merged profile is created, the number of points defining the 
profile is not verified to avoid exceeding the allocated storage for a 
maximum of 50 poiits. In the event that a merged profile should have 
more than 50 points, the overflow would incorrectly redefine neces- 
sary program locations (e.g. the near-surface temperatures, TM(1), 
TM(2), ..., would be redefined by depth values ZM(51), ZM(52), .... 
respectively, etc.) 


When the bathymetry is below the ICAPS merged profile, the inerged 
profile is augmented by ore more point defined by the bathymetry, and 
temperature and salinity values extrapolated from the deepest (last) 
depth interval in the profile. The !CAPS program does not address the 
question "How reasonable are the extrapolated values?'. Reasonabili- 
ty of these values is dependent upon how much deeper the extrapola- 
tion extends the profile with respect to the depth interval being 
extrapolated from, and tin gradient on that same interval. 


COMMENTS — ICAPS is a merge methodology that assumes the 
bathythermograph, BT, data io be the best current information upen 
which to define a sound speed profile (SSP); and, the historical data 
extending below the BT profile is the most reasonable information 
available to merge with the BT data to form a representative up-to- 
date SSP. 


1.3.1 


1.3.2 


1.3.3 


1.3.4 


1.3.4.1] 


1.3.4.2 


1.3.5 


A shallow water case is defined when the deepest depth of the input 
BT data is above 200m. This excludes a merge with historical data -- 
e.g. an operator cannot input a shallow BT profile to merge with 
historical data for a deeper water mass. In shallow areas where BT 
data is sufficient, a shallow water case sould ve specified by another 
input parameter, rather than the deepest BT profile depth. 


For a shallow water case, the salinity is set to 35 PPT. This is 
considered valid -verywhere except in the Mediterranean Sea where 38 
PPT is more appropriate. In keeping with the philusophy of making use 
of the "best" informaticn available, historical data (salinity) should 
still be utilized whenever it's available. 


Ceding of constants is more efficient when using the form 
+X.XXME+YY, rather than +X.XX*10.*#(-YY) (see subroutines SGMTST 
and SIGMAT). 


There are inconsistencies between the ICAPS program and the refer- 
ence report, "The ICAPS Water Mass History File," by Alvan Fisher, 
Jr., May 1978. 


The prograin implementation of the synthesized temperature al- 
gorithm (see Enclosure 5), is reasonable, though the approximation to 
the synthesizing algorithm is unnecessary. 


The density (sigma-t) stability algorithm described in Apperdix A 
(reference (1)) does not correspond to the implementation. The 
document states, "... salinity inversion must coincide with a tempera- 
ture inversion ...". The code assumes that if the density at depth Z. is 
greater than the density at the depth below, 25 , (starting fren the 
bottuin of the profile), then recompute the adtnity at Z, via the 
stability expression and use this "corrected" salinity value’ with the 
temperature to recalculate the density at depth Z.. There is no 
application of the criteria to adjust salinities that’ are "... within 
temperature inversions that are more than the temperature maximum 
minus 0.25°C at the lower boundary of the inversion ..." One should 
note that the stability expression was derived for a constant ternpera~ 
ture of 10°C and salinity ranges o1 30 to 40 PPT. This algorithm to 
"correct" salinity is applied whenever the density does not monotonic- 
ally increase with depth, and over the temperature canges of the BT. 
The BT temperatures can range from approximately 4°C on up to 
27°C. This wide range of temperatures suggests the need for a 
temperature varying stability algorithm —e.g. at depth Z., adjust the 
historical salinity by the increment indicated by the change in 
temperature (merge temperature minus historical temperature) to 
stabilize the water column at depth Zi -. 


Both program versions retrieve a historical profile as described in the 
report "Description of ICAPS Envi:onmental!l Data Structure" by John 
Lever, NAVOCEA"IO TN 3700-82-79, March 1979. Tie report presents 
the file structure and the retrieval! algcrithin. The following com- 
ments may enable the report to be more useful; particularly for profile 
revisions. 


1.3.5.1 


1.3.5.2 


1.3.5.3 


1.3.5.4 


1.3.5.5 


Expand the text to specifically define the ordering of profiles (water 
masses) for a geographical location. Although the ordering is intrinsic 
to the retrieval design (refer to report figure 6), the following 
conditions be'ng stated would clarify the profile data: 


When more than one water mass represent a geographical location, let 
1<N«5 be the number of profiles, then XMINT2(i), XMAXT2(i), i=l, N 
must be defined (temperature tclerances at 200m), and 


1) XMINT2(i)<XMAXT2(i) for i=1, N 
2) XMINT2(i) = XMAXT2(i-1) for i=2, N 


must be true for continuity. The exception to this is that two adjacent 
profites can have the same tolerances, XMINT2(j) = XMINT2(j+1) and 
XMAXT2(j) = XMAAT2(j+1). For this exception, there must be a non- 
zero flag (NOGL(j)) in order to consider the j+1 profile. Tie tie 
breaker is the 290-300m gradient tolerances (XMINGL(j), XMAXGL(j)) 
anu (XMINGL(j+i), XMAXGL(j+1)). When NOGL(j) is non-zero, then 
XMINGL(j) and XMAXGL(j) must be defined. Profile j is selected for 
the BT temperature gradient less than XMAXGL(j), otherwise profile 
j+l is selected. if the gradient tolerances are given for j+1, then 
XMAXGL(j) = XMINGL(j+1) must occur for continuity. In practice, the 
tolerances at j+l are only necessary to make an output message 
meaningfu: w'ien the BY gradient at 200-300m is greater than 
XMAXGL(j+1). 


On pages 3 and 9, state that the maximum number of points for a 
historical profile is 45. 


On pages 10 and 11, the references labelled "400" should be labelied 
"80" in order to agree with the programs, as do the other labelled 
references. 


On page 3, line 2 of the second paragraph, "NOPTS(N)" should be 
"NOPTI(N)"” since NOPTS is set to the NOPTI(I) value for the selected 
profile. 


On page 3, tne iext should state that NPRF«5. 
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2.1.2 


2.2 


2.2.1 


2.2.3 


2.2.4 


2.3 


2.3.1 


SIMAS 
PROGRAM ERRORS 


Conversion ervors occur involving subroutines BT and EDITBT. Metric 
inputs are converted to English units at least twice. We understand 
that this is corrected at NUSC but not at NORDA. 


Subroutine CKBT shows an error in the algorithm for modifying the 

historical data to reflect the BT profile layer depth. The DO loop at 

FORTRAN siatement number 550 may be executed with K=N15+1, i.e. 
DO 600 J=N15+1, N15 


thus, causing erroneous results. 
PROGRAM CAUTIONS 


The generated SSP arrays are allocated for 59 points. Subroutine 
INSERT adds points to the generated profile, and subroutine ENVIRN 
extends the up to 25 revised BT profile points with historical points. 
The program does not properly check on this limit. Structurally, depth 
array element Z(51) is the location of velocity C(1), the surface 
velocity which could be changed to be the Z(51) value. Also, the test 
on the 1000 ft point insert into the historical! data (line 32, ENVIRN) is 
incorrect in that DS and VS are dimensioned for 30 each. 


Subroutine BT can input up to 25 BT temperature profile points. 
Subroutine CKBT calls subroutine INSERT to augment the BT points 
with the interpolated temperature foint at 1000 feet, whenever the 
input BT profile does not have a 1000 ft. point. Since the program 
allocates space for 25 BT points, an overflow can occur, so that the 
suriace temperature, T(1), value is replaced by the overflow depth 
value, D(26). 


Whenever the BT profile is augmented by a point at 1000 ft., 
extrapolation is performed when the profile is to be extended by a 
1000 ft. point. The SIMAS program also does not consider tie 
reasonability of the extrapolation, see paragraph 1.2.).2 above. 


in subroutines ENVIRN and FORCST, there is no test that the 
corrected bottom depth value is in the deepest (last) interval which 
ends with the input fathometer (or chart) value. Consequently, 
extrapolation may occur. 


COMMENTS -- The methodology is to shift the BT SSP so that at 1000 
feet the shifted BT SSP has the same velocity value as does the 
historical prof:ie. The shifted BT SSP is then extended by those 
hisicrical prctile points be‘ow the deepest BY SSP depth. 


SIMAS assumes that the historical profile will always extend below the 
BT profile. This assumption can only be valid 1i, over each geographi- 
cal area, the associated historical profile extends down to the deepest 
bathymetry in the area, and any deepest BT depth value would not be 
deeper due to measurement errors. Though the assumption will usually 
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2.3.2 


be valid, it is too strict for a general purpose production program. 
Since the program does not verify the input data to meet this 
assumption, errored profiles may result (refer to subroutine ENVIRN 
lines 0052 through 0063). 


The criteria at 1000 feet does not assure a smooth transition below 
that depth from the bottom of the shifted BT profile and the attached 
historical points. 


Assuming that the usage of the layer depth is important, a comment 
on the SIMAS layer depth algorithm is appropriate. Subroutine LAYER 
defines the layer depth to be the shallowest depth after which the 
velocity gradient first becomes negative. This algorithm requires that 
the BT data not be noisy, and it consists cf points that only define the 
essential shape of the profile. 
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ENCLOSURE 5 
ICAPS TEMPERATURE SYNTHESIZING ALGORITHM 


Enclosure (5) to 
ODSI memo dtd 4/21/30 
Locklin to Galus 
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TEMPERATURE SYNTHESIZING ALGORITHM 


ICAPS temperature synthesizing algorithm applies when the historical profile 
extends be'ow the BT profile. The purpose is to mathc.natically merge the deep 
historical temperature profile into the bottom of the BT profile. This algorithm 
is described as follows: 


Define — the historical temperature profile to be 
(Z,, Ts j=l wyN, 
JJ 


-- the BT profile to be 
(2 ° Tat” k=ly.5Npr 


-- the merged profile to be 
(Z,, Ti), i=lyNay 


0.7006 in the Mediterranean Sea 


Tee wins 0.835 everywhere else 


Initialize the merged profile so that 
(2, Ti)Zaqs "pr, i=lyeoNpy and 
Z =Zh., j=j_ye.N 
Nets} - ‘ 
where i=j-},+1, so that depth 2H; {i=1) is the first 
° 


selected depth below the BT profile at which the synthesized 
temperatures are to begin. 


(NyeNattNy ot) 


Let T= Ty “Thy where Ty is the interpolated historical 
BT 


temperature at the bottom of the BT (Zi ). 
BT 


The synthesizing algoritl.m is given by 


(2, -Z,, )/100 
Terai hy “Ney |e 
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Clearly, the first synthesized temperature, TNBT+1 becomes TAT Gag 


as Z), approaches (merges with) the last BT depth. 


Jo 


A. Report number NOO RP-19, on page 1, presents, for the synthesized 
temperature, the equation 


-TH.+K. (K. 


NB +i 


27h )/100 


with Kaic! and Ki=A j-1 


B. The ICAPS programs use the approximating equation 


(r=} 
Typteiz 4y* Lae AT 
_(Z,-Z, 4100 
with KieA hi hey 
for i=], kK. = A@h “2p, .)/100 
Jy le leet 
oe er a azh 72h. 100] | (2, -Z,, )/100 
igtl Jo Jott lg Jo Ign! 
~ ples, Zp, 1/300 
jot! Igo! 
isi, Kauk. = A@n72p, 1/100 
eae i tg 


(Note: Z is used instead of Z , which means that the 
ey oes Nat 
oO 


product coefficient is at most off by using the 50m tolerance). 


} 
if 
| 
b 
i 
1 


ENCLOSURE 6 
TEST RESULTS 


Enclosure (6) to 
ODSI memo dtd 4/21/80 


Locklin to Galus 
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TEST RESULTS 


The Univac 1108 at NAVOCEANO was ysed to retrieve ICAPS historical data at 
two geographical locations, (42> 1'N, 9° 55' W) and (31° 3''N, 54° 36’ W), herein 
identified as cases A and B, respectively. This data served as historical data for 
execution of the CDC ICAPS and SIMAS programs. Six cases were executed by 
both programs: 


AO, BO historical data, only, defined sound speed profiles. 


Al, Bl modified historical data defined sound speed profiles. The 
modification defined the historical profile to have an average 
constant near~surface salinity. 


A2, B2 Hydrographic station data, provided by E. Hashimoto 
(NORDA 321), down to 500 meters defines a BT that is 
merged with the historical data. 


The results from these cases between ICAPS and SIMAS, though indicative, only 
represent the data for these two stations. 


1. Cases AO and BO reflect the use of Wilson's equ2tion (ICAPS) and 
Leroy's equation (SIMAS) in the generation of sound speeds. BO is the 
"worst" case and shows a difference of 0.3 ft/sec at the surface and 
becomes -0.6 ft/sec at the bottom of the profile (18144 ft.). This is 
within the relative accuracy of the models. Though both cases suggest 
that below approximately 13,000 ft the two methods may begin to 
slightly diverge. 


2. Analysis of Wilson's equation shows that, down to the 500 meters, the 
approximation of the change in velocity, V, with respect to salinity, S 
is given by 


x ~ 1.398+(3.384)10° (S-35)-(1.1)10°2T 


for T = temperature. 


Similarly, Leroy's equation yields 


x = -0.01(T -18)+1.2 


Thus, for salinities 30-40 PPT and temperatures 0-30°C 


AV ax « 1.4 m/sec for AS=|PPT 


So that, for AS = 0.1PPT, the velocity can change at most 0.14 m/sec 
(% ft/sec). 
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Cases Al and B! agree with this near-surface analysis. The largest 
salinity difference between cases AO, BO and Al, BI is at 500 meters, 
and shown in the table below. 


SIMAS 


AV 


4S 


0.6 ft/sec 


2 ft/sec 


0.14 PPT 0.5 ft/sec 


0.48 PPT 


3. Cases A2 and B2 shows the ICAPS and SIMAS merging results. Either 
method produces accoustica!ly similar profiles with the layer depths at 
the surface. The bottom (lower part) of the BT temperatures are 
within 0.1°C, so excellent agreement is expected. The following 
tables present the profile features. 


surface velocity 4958.1 ft/sec 4956.3 ft/sec] sV=1.8* 
6562.0 ft 6562.0 ft 
deep sound channel 


4922.5 ft/sec 4922.7 ft/sec 


ICAPS BT @ 984 ft = 4939.8 
SIMAS Historical @ 1000 it = 4929.5 
1.3 


surface velocity 5065.7 ft/sec 5063.2 ft/sec. | aV=2.5* 


3927.2 ft 3937.2 ft/sec. 
deep sound channel 
4896 7 ft/sec 4896.9 ft/sec. 


ICAPS BT @ 984 ft = 4991.5 
SIMAS Historical @ 1000 ft = 4990.4  4S=-l 
eas 0s i 


*The difference in aV is probably due to the SIMAS profile shift determined 
at 1000 feet. 
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